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Abstract 

We consider the problem of constructing an appropriate multivariate model for the study of the counterparty 
credit risk in credit rating migration problem. For this financial problem different multivariate Markov chain 
models were proposed. However the markovian assumption may be inappropriate for the study of the dynamic 
of credit ratings which typically show non markovian like behaviour. In this paper we develop a semi-Markov 
approach to the study of the counterparty credit risk by defining a new multivariate semi-Markov chain model. 
Methods are given for computing the transition probabilities, reliability functions and the price of a risky Credit 
Default Swap. 

keywords: counterparty credit risk; multivariate semi-Markov chains; algorithm. 

1 Introduction 

The current financial crisis has stressed the importance of the study of the correlations in the financial market. 
In this regard, the study of the risk of default of the counterparty, in any financial contract, has become crucial in 
the credit risk. For a complete treatment about credit risk we refer to the classical book of Bielecki and Rutkowski 
|l2l- Many works have been done trying to describe the counterparty risk in a Credit Default Swap (CDS) contract, 
but all these works are based on the Markovian approach to the credit risk, see for example Crepey et al. [6] or 
Ching and Ng [SJ. 

It has been shown that, the Markov chain based models, are too restrictive for the description of accurate 
rating dynamics, see for example Carty and Fons f?]. Indeed, they require that the distribution functions of the 
sojourn times in a rating class before to have a transition should be exponentially or geometrically distributed for 
continuous and discrete time models, respectively. 

In an attempt to produce more efficient credit rating models semi-Markov processes were proposed for the first 
time as applied to credit ratings in the paper by [7J; more recent results were given in ifTTl lgl. 

It is important to dispose of efficient migration models because reliable rating prediction is of interest for 
pricing rating sensitive derivatives IfTSll . IfTOl . for the valuation of portfolio of defaulting bonds, for credit risk 
management and capital allocation. 



1 



No results are available for counterparty credit risk for semi-Markov credit rating models. Such an extension 
is not straightforward as randomness in the sojourn times and memory effects are to be appropriately managed. 

In this paper, therefore, we present a novel multivariate semi-Markov model to approach the counterparty risk 
in a CDS contract. In Section 2 we investigate the behaviour of the multivariate process in the transient case 
and we derive equations for the transition probabilities with backward recurrence time dependence and reliability 
functions. An algorithm useful to perform the computations of the transition probabilities is provided. In Section 
4 we analyze the counterparty credit risk in a CDS contract. Section 5 is dedicated to a summary of this study and 
suggestions for future outlooks. 

2 Multivariate Semi-Markov Chains 

The main issue of this session is to define Multivariate Semi-Markov Chains (MVSMC), for this purpose we 
first introduce Markov Renewal chains, for more details see for example Barbu and Limnios |r|, Janssen and 
Manca {13] and Limnios and Opri^n [12j. 

Let J = (/„)„gn be a Markov Chain with values in £ = { 1 . . . . , d}, that is 

Wn+i =j\Jo = io, ■ ■ ■ Jn = i) = P{J„+i - I Jn = = Pij V G E. 

The matrix P = {pij)i.jeE is the one step transition probability matrix. Suppose that the permanence in the state is 
triggered by a renewal moment process T — (r„)„gpj, with values in N, defined by T^+i = T„ +X„ V « G N. 

X = {Xn)„fzf,} is a sequence of random variables representing the sojourn times in the n'''-state. They have a 
conditional cumulative distribution function given by 

P(:^„+i < f I ^ = i,Jn+i = j) = Fij{t) yieE. 
The couple {J, T) — 7«)nGN is said to be a Markov Renewal Chain (MRC) if 

P(7„+l = ./, T„+i -T„<k \ {Jo, To), ... , (JnJn)) = P(^+l = j,X„+i <k\Jn). 

In this case, J = {Jn)n(m is said to be the associated embedded Markov Chain. 

Further let's define the counting process N{t) = max{« G N | r„ < f }, which gives the number of transitions of 
the MRC up to time t. 

Now we are ready to introduce the Semi-Markov Chain. Let Z = (Z„)„gN be the process, with value in £ = 
{!,...,£/}, defined by 

Zk:=JN(k), k^n. (1) 

The process ([T|i is called Semi-Markov Chain (SMC) associated with the MRC (7, T), whose cumulated kernel 
is denoted by Q = {Q{k);ke N) and defined for all e E and ^ e N, by 

Q,j{k) P(7„+i = ;,X„+i <k\J„^ i). 

The element Qij{k) is the probability that the system makes next transition in state j with sojourn time less or 
equal to k given that the present state is /. 

Now, we are ready to introduce the multivariate semi-Markov chain. Let us consider a system consisting of 7 
parts, each part has values in £ = { 1 , . . . , t/}. Let's call = (/")mgn> for a = 1 , . . . , 7, the sequence of states visited 
by the a part with values in E. Denote with {T„ )„gn the sequence of transition times of the a-th component 
with state space N. Let also introduce the sequence of random variables X^ = T^^^ — T^f, for every n G N. X" is 
the sojourn time in state J". 

Let us define also the counting processes 

N"{k) = max{« G N | T,," < A:} Va = 1, . . . , 7and /t G N, 

which give us the number of transitions of the part a up to time k; we denote by N{k) = 
vector of all such numbers. 

To define the MVSMC we make two assumptions named in the following Al and A2 



{N^{k),...,Nl'{k)), the 
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Al The process J ~ (-/")a=i, ■ ,y is a multivariate Markov chain in the sense described below. 

Given the component a , we denote the vector of all components except of a by the symbol — a = (l,2,...,a — 
l , a + 1 , 7). We introduce the marginal one step transition probability for the multivariate Markov chain 
J for all i e N as 

I CT(7«,/z<A^«(.)),(T(y,-«,/z<N-«(.) + l)) 

= lP(-^^a + i = I -^"a (s) ; -^N "a (,.) ) =: Pj^^^j . j («) (2) 

where G(J",h < N"{s)) denotes the natural filtration of 7" and a{J^",h < N^"{s) + 1)) is the natural 
filtration of the multivariate J^" process. We assume also that the process is time homogeneous and then, 
the transition probabilities do not depend on time s, so we have 

A2 For every a, j3 = 1, • • • , 7, the sequences of sojourn times (^")„gn and (X,f )„gN are independent in the sense 
of formula (|4]) here below: 



F-{k;s) := P(X«„(,)^i<fc|(T(7;,^,^j,xJ,(^^),/.<.,^ = l,...,7) 

= p(x;+i<fc|7« = /„)=^;-:« v.eN, (4) 

where a{J^p^^yX^p^^yh < 5,j3 = 1, ... ,7) is the natural filtration of the multidimensional (J,X) process. 

ice) 

Then, in our model we suppose that the sequence of sojourn times {X„ )„eN depends only on the visited state 
of the a -component J". In the second step we use the independence of the a part sojourn time from the trajectory 
of the other component of the system. 

The discrete time semi-Markov kernel for each part a of the system can be formally defined as follow 

Qi jik-s) P(y"a(^,)_|_i ^ j,X^a^,)^i < k \ 7^,^^,^ = ; 1 , . . . , 7^^^^^.^ = ly), (5) 

for all i G N, and where i = (i'l , . . . , iy) is a vector in E"^ and j is an element of E. By hypothesis we have that the 
semi-Markov kernel does not depend on time s, that is 

Q?jik;s) = Q^jik) ysen. (6) 

Notice that the semi-Markov kernel introduced in formula (|5]l can be written as 

Q?j{k)=p?jFrjk), (7) 

where it is implied that i ~ Jn(.s)+i and j — Jn''{s)+i, for all i e N. 

In what follow, it can be useful to speak in terms of probability to have a transition exactly at a certain time, 
hence we define 

^i"; = ^ i-^N" (.«) + !=■/' ^N" (i) + 1 = ^ I -^W 1 (s) = ' 1 ' ■ ■ ■ ' ■^iVr(i) ^ ' r) 

= e?.(fc+l)-Gr/fc)-:g?.W, (8) 
for all i e N. Notice that the transition probabilities (j2ji can be expressed in terms of (^fy) as 

prj=i^m (9) 

k=0 
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Figure 1: The trajectory of the double-component system is shown as a function of time. In the picture sojourn 
times, transition times and backward recurrence times are shown. 



We define, for each component, for every state ; G E and time e N, the unconditional sojourn time cdf: 

:= P(X«+i <k\J^ = i) = P(X«„(.^.)+i < k I = /„) 

= V.eN, (10) 

We notice that in this case we have //" = F", that is 

= (11) 

The evolution of the multivariate (J,X) process can be described as follow, given an initial state the next state 
occupied by the system is determined according to the evolution of the multivariate Markov chain while the sojourn 
time in the present state, of every part of the system, is determined according to the joint distribution of X. In the 
bivariate case an example of trajectory is shown in Figure ([TJ. 

Now, we can introduce the MVSMC associated to the kernel {Q")a=i.- - .y- We denote the MVSMC by 
Z = (Z")a=i,...,y where Va = 1, . . . , 7 

Z«(^):==y«„(^) ykeN. (12) 
The transition function of the MSMC for a-th component is defined by 

4>?.(^)=P(Z«(^)=;|Z(0)=i), (13) 

for all ieE^JeE and keN. 

In the semi-Markov environment, the Markovian property is preserved only at transition times, in general we 
have to take into account for the history of the process. Then transition probabilities depend on the elapsed time in 
a given state, see for example D'Amico et al. lIS). 

In general when we look for transition functions of semi-Markov model we need to introduce the concept of 
backward recurrence time. 
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The backward recurrence time for any component a, denoted by B", specifies at any time the age of the state 
of the a-th part, that is 

B«:=f-r^«(,) Va = l,...,7. 
In other words, Bf gives the time from the last transition of a-th component. 

The one step transition probabilities for component a, by tacking into account for bacward recurrence time, is 
defined by 



= i'^iv«(.v)+i =k + va\ Jn(.v) = i , Tn(,) =s-\ , Tn(,)+i > s) qij{va,k) 

1 



where, Va stand for the backward of component a at time s, i = (/i , . . . and Jj^^.j = {jj^ . ^) are vectors 



in £1', V = (vi , . . . , vy) and Tn(.) = (r^^ ( J , . . . , T^J^^.j ) are vectors in W. 



Proposition 1. The one step transition probabilities for component a are given by 

a, /■.^(A: + v„)-^.«(fc + v„-l) 

Proof. By applying Bayes rules we get 

+ l = ^''^"«(.s) + l = ^ + >-'0! I Jn(.s) = i J TN(i.) = i - V , TN(i)+i > 
= P(^"a(i) + 1 = ^+ I /"a(,,)_|_i = , Jn(i) = i , Tn(,.) =S- \ , Tn(,.)+1 > s) 

X = i I Jn(.s) = i > Tn(.,) = i - V , Tn(,.)+i > i). (15) 

Then, by using assumptions Al and A2 we obtain 

= j^^m(s)+\ = ^ + '-'a I Jn(.s) = i : Tn(s) = i - V , T]M(j)+i > i) 
= IP(^iV« (.() + ! = ^ + >-'0! I ■^"a(i) = 'a , r^a(i) = - t-'a , ^N(i)+1 > 

X JNM=i), (16) 



and by using the definitions and formula ( 11 1 we get 



= ■U^m{s)+\ =k + ^a\ Jn(.s) i , Tn(s) = 5 - V , Tn(s)+i > 
^■:(fc + va)-^:(fc + va-l) ^ f;-:(^ + v«)-^-:(^ + Va-l) 



□ 



Notice that if the backward process is zero [s) — {va — 0) we recover 



Formula (14i reveals that assumption Al and A2 imply that the one-step probabilities q'^ are affected by the 
duration in the state only of the component a. That is, the backward values of the other component does not affect 
the kernel q". 

Anyway, as we will show, the backward values of each component affect the transition probabilities of each 
other components. 

Now we want to discuss the evolution equation for the MVSMC with backward recurrence times. First of all, 
we define the transition probability for the component a as 

4>?.^(v,Ma,^):=P(Z«(^)-j„,B«(^)-M„ |Z(0)=i, B(0)=v). (18) 

The following result consists in a recursive system of equations for computing the transition probability func- 
tions in a bivariate system. 
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Proposition 2. Suppose that the system is composed of two components, i.e. 7=2, and such that hypotheses Al, 
A2 hold true. Then, for a// i, j G E'^, v, u G and G N, we have 

l-H}{k + vi) 

k r+V2 
T=l/i,/2G£ H'=0 



(19) 



(V2) 



+ L L L *(/,.h)J,((>^'0),M2,fc-T)*[,.^.,.^)_,^((vi,V2),W,T)4. _,.^).,^(V2,^ 

t=lli,l2eE w=0 

Proof. We show the resuh for transition function of component 1, all what follows holds for component 2 by 
symmetric arguments. 

For all i,j G E^, v,u G and A: G N, we make a partition of the sample space in two parts; no transition up to 
time k and at least one transition; therefore we have 

4,. {y.uuk) = V{z\k) = h,B\k) = Ml I Z(0) = i , B(0) = v) 



= P(zi(^)=ji,Z?'(^) = Mi,ri' >fe|Z(0)=i, B(0)=v) 

k 

+ W = juB' {k) = u,ji=x\ Z(0) = i , B(0) = v). 

T=l 

Let us consider the first term on the right hand side, using Bayes rules we obtain 

P(Zi(A:)=yi,Bi(A:) = Mi,r/ >A:|Z(0)=i, B(0)=v) 
= n2\k)=juB\k) = u, I r/ Z(0)=i, B(0)=v) 
X P(r/ > k I Z(0) = i , B(0) = v) 
l-//i(fe + vi) 



'''■lil^{«I=*:+vi}" 



i-H,;(vi) 



Now, let us consider the second term on the right hand side of equation pO| , we have 



' (fe) = ;i (fc) = «i , 7^1 = T I Z(0) = i , B(0) = v) 

= L P(z'(^)=ji,B'(fe) = «i,r/=T,Z(T)=l, 

h.heE M'=o 

B2(T)^w|Z(0)=i, B(0)=v), 
where = t imply Z?' (t) = 0. Then by applying Bayes rules we obtain 

P(Z1 {k) = juB' {k)^ui,Tl^T\ Z(0) = i , B(0) = v) 

= £ £P(z'W = ;i,B'W = «i |z(T)=i,r/ = T,B2(T). 

ll,l2<^E w=Q 

X P(z1(t) = /i, = t I Z(0) = i,BHO) = vi,b2(t) ^ w,z\t) = h) 
X P(z2(t) = /2,b2(t) = w I Z(0) = i,B(0) = v). 

Here, by using the time homogeneity, the first term in the right hand side become 

P(Z1 {k) = y 1 (^) = Ml |Z(t) = 1, r/ = T, b2(t) = w) = <J>[,„,^) ((0, w) , mj , ^ - t) , 



w) 



(20) 



(21) 



(22) 



(23) 



(24) 



the last term, by definition, is 

P(Z2(T) = Z2,B'(t) = W I Z(0) -i,B(0) =V) =4>2^ ,^) ,^((V1,V2),W,T). (25) 

Finally, let us consider the second term, we have 

P(z1(t) = Zi,r/ = t I Z(0) =i,Bi(0) = vi,b2(t) = w,z2(t) = Zz) 
= P(y^.(,) = = T I Z(0) - i,Bi(0) - yi,B2(T) = w,z2(t) = /2) 

X P(4,(^)-Zi |Z(0)-i,B'(0)-vi,B2(T)^vv,z2(T)-/2). (26) 
Here, by using assumption A2, we obtain 

^(T^Hr) = ^ I = ^i'Z(O) = i,B'(0) - vi,b2(t) = w,z2(t) = ^2) 

= P(r;.(^) = T I = /i,Z' (0) = /i,b'(0) = vi) 

= iP(V(r) = ^|2 (o) = 'i''B (o)-vi) = i:^F^) 



furthermore, by using Al, we get 



Then 



P(zi (t) = IuTI=z\ Z(0) = i,B\Q) = vuB^it) = w,z\x) = h 

j^.;(T + v«)-j^. ;(T + vi-l) 
1 



By substituting the formulas ( |24|25|29] l in ( 23 1 we get 



(27) 



P(4,(^) - /i I Z(0) = i,Bi(0) = vy,B\x) - vv,z2(t) = Z2) 

P(4'(r)=^i I JN(0)=i)=W./r (28) 



Pu, ='?i./,(v'i,T). (29) 



' (^) = ;i (^) = «i , 7^1 - T I Z(0) - i , B(0) - v) 

T+V2 
hdl^E vv=0 

Finally, by insert this last equation in pOb we get the result. □ 



2.1 The algorithm 



To apply the model we must solve the system of equations (19). The main steps of the algorithm are here 
described. As first step, by putting A; = 1 we can have only two possible values for the final backward process, that 
is Ml = or Ml = vi + 1. 

Let consider first the case k—\ and mi = for the first component. 

By considering that mi = implies that l.^ui=\+vi} = we can rewrite the first equation of system (19) as 
follows: 

O?,, (v,0,l)= L^(',.,,)./,(^i-lK/,.,)./,((^i.^2),vv,l)<I>|,^y_. ^ (31) 
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Now, notice that 4>|^^ ^.^ ((0, w),0,0) — 5ij^ and by substitution we get 

1+V2 

^L.(^'0,i)= 4„„),,,(vi,i)5,„. L«>^,-„„),,((vi,v2),w,i) 

/lG£ 

The second case to deah with is A; = 1 and mi = 1 + vi for the first component. 
In this evenience we have simply that 



, 1 -//,.! (1+vi) 

4>^(v,l+vi,l)^g,,. • (33) 



In fact the second addend on the right hand side of the first equation in ( 19 1 is zero because it contains the term 
'2)71 ((^'^)' ^ + ^iiO) =0- The reason of this last equality is that, due to the fact that k = Q, the initial and 
final backward values of the first component should be equal, but here we have ^ 1 + vi because vi cannot be 
negative. 

This first step has been executed for the first component, a symmetric argument gives similar results for the 
second components. 

Observe that at time k = \ the transition probability of one component is affected by the state of the other 
component but not by the duration in the state of the latter 

As second step, by putting ^ = 2 in ( 19 ) we have the following equation: 

2 T+V2 

+ L E E '?(V.2)Ji(^i'^)'^L'2)i2((^i'^'2)'^'^H,A),;i((0,vv),Mi,2--T). 

T=l/i./2G£n'=0 

When k ~2 and the initial backward of the first component is (0) = vi , then B' (2) e {0, 1 , 2 + vi }, conse- 
quently we have three cases to analyse: (k = 2 and mi = 2 + vi, k = 2 and u\ — \,k — 2 and u\ — 0). 
If A; = 2 and ui ~2 + vi we have 

I+V2 

+ E E4^.2)yl(^l'l)*(.^<2)i2((^''^2),w,l)4>|,,_,,).^. ((0,w),2 + vi,l) (35) 

/l,/2G£ H'=0 
2+V2 

+ E E 4..2),J^i'2K-,.2)./,((^'i'^2),vv,2)O(^,,)..^((0,vv),2 + v,,0). 

Zl ,/2G£ H'=0 



Note that the third added on the rh.s. of (|35|) vanishes due to the fact that 

('l,'2)Jl' 



4>}, . ((0, w), 2 + yi,0)=0 being 2 + VI 7^0. Then, It result that 



1+V2 ' (36) 

+ E E 4,i2),/,(^l'lK^Z2)J2((^'l'^'2)'^'l)<,/2).;■.((0,vv),2 + Vl,l), 

/l,/2G£ M'=0 

moreover note that the the third factor on second addend on the right hand side vanishes, in fact by definition 

'^{/i,/2u((0,>v),2 + vi,l) =P(Z1(1) =2 + yi | Z(0) =1 , B(0) = (0,vv)) =0 (37) 
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and finally we get 

l-///(2 + vi) 

^ , (V, 2 + VI , 2) = , (38) 



as we can expect, the probability to go in 2 unit times from an initial backward state vi to a final backward state 
vi + 2 is the probability to make no transition in the 2 unit times. 
The second case is when k ^2 and mi = 1 where we have 

4>K (v, 1,2) = g/i,il{i=2+vi} ';'^r"^ ^ 

1+V2 



' L L^('^,2)7:^'l'2)4>f,-,,„),z,((vi,V2),vv,2)4>J,,_,,)..^((0,vv),l,0). 



(39) 



This last equation, due to the fact that I{i=2+vi} — and 



^(h h) u ((0>>v)> 1-0) = 0, simplifies in 



(V,l,2) = ^ 'f ),^(vi,l)02. ,^),^((V,,V2),W,1)0|,^_,^).^ 



now, form Step 1 (cf Eq.( 33 1), the third factor on right hand side is given by 

'^Ih .h).h ( (0, vf) , 1 , 1 ) = 5,, ( 1 - //,; ( 1 ) ) = 5/, ( 1 - //,\ ( 1 ) ) , (40) 
then, by substitution, we finally get 

(v,l,2)=^[„,,^) .,(vi,l)(l-//j,(l)), (41) 

as we can expect, the probability to go in 2 unit times from an initial backward state v\ to a final backward state 1 
is just the probability to make one transition on the first unit time and then wait one uint time. 
The last case is A: = 2 and u\=Q where we have 

<5i^,(v,0,2)^5,,.I{o.2+va^^ST^ + 

+ L L^(V.2).'l(^'''l)*(M,-2)./2((^l'^2),VV,l)<I>(V,,)..^((0,Vv),0,l) (42) 

lihf^E w=0 



Remark that I{o=2+vi} — ™d <I>|^^ ((0, w),0,0) = make (42 1 equal to 

1+V2 

OJ,. (v,0,2)=0+ £ ^ ^[,^..)_,^(vi,l)4>2. _,,),,^((vi,V2),vv,l)Oj,,,y,^.^((0,w) 
hh'iE vf=0 

2+V2 

,!2),'7((^i'^2),vi',2)5(/j 

/1./2G-E »-=o 

1+V2 

= E E4^,2yl(^l'l)*k<2)./2((^l'^2)'^'lK^y,yl((0,vv),0,l) (43) 

/l,/2G£ M.'=0 

2+t'2 

+4i,'2),.i(^i'2) ^ ^ 4>f,__,^)_,^((vi,V2),w,2) 
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which are all known terms from Step 1 . A symmetric argument gives similar results for the second components. 
At this point, proceeding this forward algorithm we can get O"^^ {\,Ua,k) by knowing <I>?^^ {\,up,s) for all i < A; 
andj3e{l,2}. 



3 The 2-component reliability model 

Let us consider two partitions of the state space E of the two components as follows: 

E = U'[jD\u'r]D^ =(0; 
E = U^\jD\u^f]D^=(I). ' 

The subset (U^) contains all good states in which the component 1 (2) is regarded as well working and the 
subset (D^) contains all the bad states in which the first (second) component is not well performing. 
Let us formulate an additional assumption named A3. 

A3 All states in (D^) are absorbing for the component 1 (2). 

Notice that, the assumption A3 can be relaxed with easiness and it is possible to execute the following computation 
in the general case in which the first (second) component alternates between the set t/' and D' (JJ^ and D^). 
Anyway we adopt A3 because the presentation is easier and furthermore the application we discuss in next section 
falls well within this case. 

One of the most useful indicators is the reliability function. Here we define the reliability of the system as 
follows; 

Ri,i^{{vuvi),k) := V{Z\k) e U\Z^{k) e f/2 I Z(0) =i,B(0) = v). (45) 
In what follows we need also of the marginal reliability functions for a single component defined as: 



4.,2((^i>^2),^):=P(Z'WGf/MZ(0)=i,B(0)=v) 

<,„((vi,V2),^) :=P(Z2(^) e f/2 I Z(0) =i,B(0) = v) ^ 



The following results expresses the link between the marginal reliabilities and the system reliability in our 
model: 

Proposition 3. Suppose that the system is composed of two components, i.e. J —2, and such that hypotheses Al, 
A2 and A3 hold true. Then, for all i, j G -E^, v, G and A: G N, we have 



Ri,i,{{vuV2),k}^RliA{vuV2),k)RlUi^uV2),k) (47) 



where for a G {1,2} we have 



Proof. By using Bayes rule we get 

7?„,,((vi,vi),fc) = p(zi(fc) eu\zHk)eU^\ z(0) =i,B(0) = v) 

= ¥{z\k) G I Z^{k) G C/^Z(0) = i,B(0) = y)V{Z^{k) G \ Z(0) = i,B(0) = v), ^ ' 

and by using formula ( [T8] l, the definitions of the semi-Markov and backward processes and assumption A3, we get 
in 



By using assumption Al and A2 we have 



(50) 



J2eu^ ' (51) 

= L L <PljJy,;k)=Rl{{vuV2),k)Rl{{vuV2),k). 
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□ 



Remark 1. The precedent result allow us to express all the reliabilities in terms of the transition probabilities. 

Remark 2. If we consider a model where the two components are supposed to be independent, then, the reliability 
of the single component would be modelled with a standard univariate semi-Markov chain and we would get 



If the two components are supposed to be independent, the product of these two reliabilities should be equal to the 



by one is an indication of the correlation between the two components. 

4 Counterparty Credit Risk in a CDS contract 

In the financial market all subjects are exposed to the default risk. Then, in any financial contract we have to 
take into account for the risk of default of the our counterpart. Counterpary credit risk is in general 'the risk that 
a counterpart of a financial contract will default prior to the expiration of the contract and will not make all the 
payments required by the contract' (cf. Pykhtin and Zhu [14]). 

We would like to study the counterparty credit risk in a Credit Default Swap (CDS) contract. In this work we 
would like to emphasize the difference between the CDS contract with and without consider counterparty risk, we 
will call these two cases risky CDS and risk free CDS respectively (Crepey et al. |6 1). 

Let us consider a firm C, supposed to be defaultable, emitting an obligation (or bond) on one money unit at the 
time with maturity time T. Let us also consider a bondholder A (or protection buyer) supposed to be risk free in 
all what follows. The possible financial scenarios are 

• If C has not been default until T, it is able to pay the money due to the bondholder A. 

• In case of C's default before or at the maturity date T, it will be able only to pay a fraction (recovery rate 
Pc) of the face value of the obligation to A. 

For these reasons the bondholder A is looking for protection against the loss in case of C's default. Let us consider 
a third financial subject that we will call generically as protection seller B. A risk free CDS is a contract which 
obligates A (protection buyer) to pay a fee to B (protection seller supposed to be risk free) in change of protection 
against the default of the reference credit firm C. The cash flows of a risk free CDS are 

• A pays to B a stream of premia with spread k, from the initial date until the occurence of default event or the 
maturity date T. 

• In case of default of C, B has to cover the loss of A. Then B has to pay 1 — pc unit of money to A. 

The value of the spread is evaluated in order to guarantee, that the contract has value zero at the inception time. 
We assume that the payment of B to A is made at the same time of the default event. 

Let Tc be the time of default for the credit reference firm C. From the above discussion, we can directly write an 
expression for the cash flows and price process of the risk free CDS contract. The In Cash Flows process from the 
perspective of the bondholder A in the risk free CDS is given by 



Rlivuk) :=P(Zi(A:) e \ Z\0) = iuB'{0) = vi) 

4(v2,fe) :=P(Z2(^) e f/2 I Z2(0) = /2,B2(0) ^V2). 



(52) 




^Tc(l-Pc)I{r<Tc<n' 



(53) 
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where j3 is a discount factor. The Out Cash Flows process is given by 

T 

-^:^A.I{.v<rc}- (54) 

s=t 

Then, the discounted value of the risk free CDS with maturity T is 

T 

Prprit) = -KY^Iisl{s<rc}+^A^-pc)l{,<rc<T} 

.s=r 

and its price process is given by P, = E([/:i7'(f)]. The subscript t, here and after, indicates the information that at 
time t the process is still in one of the Up states. 

Remark 3. The price of a risk free CDS can be evaluated with a single component reliability model where the 
rating of the only defaultable subject is modeled via a standard univariate semi-Markov process (see for example 
D'Amico et al. [W]). 

A risky CDS is a contract which obligates A (protection buyer) to pay a fee to B (defaultable protection seller) 
in change of protection against the default of the reference credit firm C. The cash flows of a risky CDS are 

• A pays to B a stream of premia with spread k, from the initial date until the occurence of default event or the 
maturity date T. 

• In case of default of C, if B has not defaulted, B has to cover the loss of A. Then B has to pay 1 — pc unit of 
money to A. 

• In case of default of B, if C has not defaulted, the contract is stopped with a Close-Out Cash Flow (cf Crepey 
et al. fE\). In this work we assume that the two parties according on a termination of the contract with a 
terminal cash flow paid to A, positive or negative, depending on the value of the risk free CDS computed at 
the time of default (cf. Brigo et al. |3 1). 

• If B defaults at the same time as the firm C, B will be only able to pay to A a fraction (recovery rate pB)of 
the loss of A, namely Pb(1 — pc) unit of money. 

The value of the spread is evaluated in order to guarantee, that the contract has value zero at the inception time. 
We assume that the payment of B to A is made at the same time of the default event(s). The possible loss of A for 
the joint default event is an effect due to the counterparty risk. 

Let us introduce Tb, the time of default for the protection seller B. The In Cash Flows process for the risky 
CDS is given by 

Prci^ - Pc)\<tc<T}[\tc<Tb} + Pb\tc=Ts}]+ PTBl{,<rs<{TArc)}PBPtB (55) 

here j3 is a discount factor and with P.+ we denote the positive part of the price process for the risk free CDS. The 
Out Cash Flows process is given by 

T 

^ ^Y.Ps^{s<(tcATb)} - l^rB^{t<TB<{TMc)}^TB (56) 

s=t 

where Pr stands for the negative part of the price process for the risk free CDS. Then, the discounted value of the 
risky CDS with maturity T is 

T 

PtTtrit) = -Kj^lisl{s<(rcATB)}+Prc{^-pc)\<Tc<T}[MTc<TB} + 

+ PBl{rc = rB}]+PTB^<TB<{TA.Tc)}ipBPtB~PtB)^ (57) 

and the price process for the risky CDS is 11, = EJ;r7-(f)]. 
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4.1 Pricing a risky CDS and CVA evaluation 

In this section we are going to apply the 2-component reliabihty model to the study of the counterpart risk in 
a CDS contract. In particular our goal is to price a risky CDS contract and to derive an expression for the credit 
value adjustment (CVA) which can be seen as a measure of the counterparty credit risk. 

In order to price a risky CDS we should be able to evaluate the expected value of the indicator functions in ( [57] l. 
The following result concern the evaluation of the expectation. 

Proposition 4. The price of a risky CDS under the natural probability measure is 



mT{t)^~KY^^ E Pr(T = /l)+ f; f j3,,^(l-pc)P,(Tc = /lc,Tfi=M 

s=t h=s+\ hB=t+\hc=t+l 

T 

+ E j3/,c(l-Pc)PslPr(Tc = /Jc,TB=/zc), 

hc=t+l 



(58) 



where T = Tc A Xb, 



and 



,(T = /i)= £ V,{Xc^h,XB=hB)+ E V,{Xc^hc,TB^h), (59) 

hB=h hc=^^ 



P, (tc - /^c, - - E E (^i.y i^'c - 1 ) - <v (he) ) «v (hB - I) - R^M) . (60) 

Proof. The result is a direct consequence of formula ( [57] i but it remain to prove that 

P, (tc = /^c, E E (-^Lv C'c -l)-Rt{hc)) {hB -l)-Rt{hB)). (61) 

We will show the result for f = for the sake of simplicity, the general case is a direct consequence. First of all we 
notice that 

Po(tc = /ic,Tb = M- E E IP(2^(M=AZ''(M =-0 12(0) =i,B(0)=v). (62) 

iGC/2 vGN2 



(63) 



Using the properties of the probability measures, we get 

P(zC(/2c) D,Z^(/!b) =. D I Z(0) - i,B(0) = v) = 1 - [P{Z^{hc) € U \ Z(0) = i,B(0) = v) 
+P(Z«(/ib) e U I Z(0) = i,B(0) - v) -P(zC(/;c) e C/,Z^(/zb) G U \ Z(0) - i,B(0) - v)] . 
The last expression can be rewritten in terms of reliabilities as 

P(zC(/ic) = £>,Z^(/7b) = D I Z(0) = i,B(0) = v) = {Rf^^{hc~~l)-Rf^,{hc)){R?,ihB-l) -Rl,{hB)), (64) 
which concludes the proof. □ 

Remark 4. The precedent result gives an expression of the price of a risky CDS as a function of the reliabilities, 
that is the transition probabilities for the bivariate semi-Markov chain (cf. remark^. 

Remark 5. The difference between the price of a risky CDS and the price of a risk free CDS is has a particular 
importance, indeed it is a measure of the loss of value a CDS contract undergoes due to the counterpart credit risk. 
This difference is called credit value adjustment ( CVA). The credit value adjustment process (CVAt) is defined by 

CVA, =P,-Yl, fort<ZB, (65) 

it measures the loss of value of the CDS contract. We notice that, in our model, the CVA process can be totally 
expressed in term of the reliabilities. 
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5 Conclusions 



This paper proposes a multivariate semi-Markov chain model in discrete time. The multicomponent system is 
analysed in the transient case by giving methods for computing the transition probabilities and reliability functions. 
The numerical solution of the equations is made available by means of a recursive algorithm. 

The results are applied to the evaluation of the risky credit default swap contracts and they allow the attainment 
of an explicit formula for the price of a risky CDS and for the credit value adjustment process. 

Possible avenues for future developments of our model include: 

a) application to real data on credit rating dynamics; 

b) asymptotic properties of the multivariate semi-Markov model; 

c) construction of a multivariate reward model for the credit spread computation. 
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